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Recent work suggests nontrivial relations between generalized parton distributions on the one 
hand and (naive time-reversal odd) transverse momentum dependent distributions on the other. 
Here we review the present knowledge on such type of relations. Moreover, as far as spectator 
model calculations are concerned, the existing results are considerably extended. While various 
relations between the two types of parton distributions can be found in the framework of spectator 
models, so far no nontrivial model-independent relations have been established. 

PACS numbers: 12.38.Bx, 12.40.-y, 13.85.Ni, 13.85. Qk 



I. INTRODUCTION 

For decades the partonic structure of the nucleon was 
almost entirely discussed in terms of unpolarizcd and po- 
larized forward parton densities, which merely depend on 
the longitudinal momentum of the respective parton. A 
much more comprehensive picture of the nucleon struc- 
ture, however, can be obtained by considering two other, 
more general types of parton distributions: first, general- 
ized parton distributions (GPDs) entering the QCD de- 
scription of hard exclusive reactions on the nucleon (see, 
e.g., Refs. [ILIIIIII]); second, transverse momentum depen- 
dent parton distributions (TMDs) entering the descrip- 
tion of various hard semi-inclusive reactions (see, e.g., 
Ref. 0,11,1). Not only a large body of theoretical work 
on these types of parton correlators appeared during the 
last decade, but also new high luminosity particle ac- 
celerators nowadays allow one to explore such intriguing 
though complicated objects experimentally. 

On the TMD side the so-called naive time-reversal 
odd (T-odd) parton distributions are of particular im- 
portance, because these objects can give rise to single 
spin asymmetries (SSAs). Single spin phenomena were 
measured in hadron-hadron collisions at FermiLab 0, [1] 
and at RHIC [9|, [l0( , as well as in lepton-hadron collisions 
by the COMPASS ConaborationJll|Jll, the HERMES 
Collaboration [H EI EE El ElSiEl, and at Jef- 
fersonLab [2(j. In general, the theoretical description 
of such observables in the framework of QCD has been 
and still is a challenge for the theory. However, a cru- 
cial step forward was the observation that time-reversal 
invariance of the strong interaction does not forbid the 
existence of T-odd TMDs [U At least for the SSAs 
in lepton-induced reactions a QCD-factorization formula 
containing T-odd TMDs has been put forward in the 
meantime [H, [H, HU . 

An important object in this context is the T-odd Sivers 
function [H [2?J , denoted by in the nomenclature of 
Ref. (28|. The Sivers function quantifies the SSA related 
to the transverse momentum dependent distribution of 
unpolarized partons inside a transversely polarized tar- 



get. Experimentally can be studied, e.g., by measur- 
ing a transverse SSA in semi-inclusive deep- inelastic scat- 
tering (SIDIS). Extractions of the Sivers function from 
existing SIDIS data E1IME1E1 
performed [U, HI Ullla IS | . 



have already been 



An intuitive picture of various transverse SSAs, in 
particular, also of the Sivers effect, was proposed in 
Refs. [H, [36J. That work suggested for the first time 
that there may be a close connection between a certain 
GPD, typically denoted by E (see, e.g., Ref. Q), and the 
Sivers asymmetry. In order to generate the Sivers effect 
in SIDIS two ingredients are required in the picture of 
Refs. [35l |3^| : first, the impact parameter distribution of 
unpolarized quarks in a transversely polarized target has 
to be distorted. Such a distortion is directly connected 
to the GPD E. Second, the fragmenting struck quark in 
SIDIS has to experience a final state interaction with the 
target spectators. 

Though the work in [35|, Hi| provided a very attractive 
explanation of the origin of the Sivers effect, it did not 
provide a quantitative connection between the GPD E 
and fi T . For a particular moment of the Sivers function 
such a connection was later on established in a perturba- 
tive low order calculation using a simple diquark specta- 
tor model of the nucleon [37j | . Another step forward was 
made by comparing the correlator for chiral-odd GPDs in 
impact parameter space on the one hand with the correla- 
tor for chiral-odd TMDs on the other [H]. The work [H| 
suggested, in particular, a relation between the second 
leading twist T-odd quark TMD, the Boer-Mulders func- 
tion [28j], and a certain linear combination of GPDs. 
This possible connection between and GPDs was af- 
terwards discussed in more detail in Ref. [3{|. In the 
framework of the diquark spectator model quite recently 
another relation between the GPD E and the Sivers func- 
tion was obtained [4oT ] . which is similar to the one found 
in |22|. 

Despite these developments so far no nontrivial model- 
independent relations between GPDs and TMDs have 
been established. (An attempt in this direction was 
made, e.g., in Ref. 36].) If one considers for instance 
the Sivers function, the crucial problem lies in the fact 
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that there is apparently no model-independent factoriza- 
tion of fyj, into a distortion effect, described by the GPD 
E in impact parameter space, times a final state interac- 
tion. 

The manuscript is organized in the following way. In 
Sec. |TI] we give our definitions of the GPDs, both in 
momentum and in impact parameter space, and of the 
TMDs. In particular, we also include the gluon sector 
which so far in the literature has not been discussed at 
all in the context of possible relations between GPDs and 
TMDs. Section HTT1 is devoted to model-independent con- 
siderations of the relations between GPDs and TMDs. 
Here the current knowledge on this point is summarized 
and the main difficulties are presented. Also new possi- 
ble relations for gluon parton distributions are provided, 
which later on in the manuscript are investigated in a 
model calculation. Section HVl describes model results on 
relations between GPDs and TMDs, where we exploit 
two perturbative models for the target: the scalar di- 
quark spectator model of the nucleon, and a quark target 
model treated in perturbative QCD. The latter, in par- 
ticular, allows one to study the gluon sector. The model 
results of the various parton distributions, calculated to 
lowest nontrivial order in perturbation theory, show rela- 
tions in accordance with the considerations in Sec. IIIIL In 
the framework of the two models also on the quark sec- 
tor new relations are established. Some of these relations 
contain the results of Ref. [13] and of Ref. [4(| as limit- 
ing cases. In this section we also argue that even in the 
context of spectator models certain nontrivial relations 
between GPDs and TMDs are far from being obvious if 
one considers higher orders in perturbation theory. Our 
summary is given in Sec. |VJ The results on the var- 
ious parton distributions in the two target models are 
collected in two appendices. 



II. DEFINITION OF PARTON DISTRIBUTIONS 
A. Generalized parton distributions (GPDs) 

We start by presenting the definitions of the GPDs. 
Unless stated otherwise we follow here the conventions 
of Ref. [2| . The momenta of the incoming and outgoing 
nucleon are given by (see also Fig. [1} 



p = P-±A, p' = P + ±A. 



(1) 



and satisfy p 2 = p' 2 = M 2 , with M denoting the nucleon 
mass. The GPDs depend on the three variables 



A H 



t = A 



' p+ > - 2 P^ 
where the light-cone coordinates are defined by 



(2) 



v ± = ^(v°±v 3 ), v T = (v\v 2 ) (3) 

for a generic 4- vector v. In a physical process the so- 
called skewness £ and the momentum transfer t to the 




FIG. 1: Kinematics for GPDs. 



P 



nucleon are fixed by the external kinematics, whereas x 
is typically an integration variable. It is convenient to 
define the following tensors, 



°T — Hi 



Jj _ , — Hi 



Moreover, we use the conventions 



£ 0123 = 1 i 



■ 1 2 3 

75 = »7 7 7 7 > 



t^ v = U^^f - 7^7 M ) 



(4) 

(5) 
(6) 
(7) 



The quark GPDs are defined through the light-cone 
correlation function 



*W{-\z;\z)l>{\z) |p;A) 



z+=q+ ' 



(8) 



where A and A' respectively characterize the helicity of 
the nucleon in the initial and final state. The object V 
is a generic matrix in Dirac space. In ([5]) a summation 
over the color of the quark fields is understood. Fur- 
thermore, in a perturbative calculation of the correlation 
function only connected graphs have to be taken into ac- 
count. The correlation function in ((5|) also depends on 
a renormalization scale which we disregard throughout 
this work. The Wilson line W in ([8]) , connecting the two 
quark fields and ensuring color gauge invariance of the 
correlator, is given by 



W{-\z-\z) 



z + =0+ 
zt— Ot 



[0+,-iz-,0 T ;0 + ,iz-,0 T ] 



2 ,UJ '" '2 
/ HI/ 2 )*" \ 

= Vexp(-ig dy-t a Ai(0+,y-,d T )), (9) 

V J-(l/2)z- J 

with V denoting path-ordering and t a representing the 
Gell-Mann matrices. For three particular matrices T in 
© one obtains the leading twist (twist-2) GPDs: 

F q (x, A; A, A') = (x, A; A, A') 



2M 



E«(x,S,t))u(p, A) 



(10) 
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F"(x, A; A, A') = F^ + ^ (x, A; A, A') 

+ ^T&(x,t,t))u(p,\), 
F«< j (x, A; A, A') = F q[ia3+ ^ ] {x, A; A, A') 



(11) 



2M 

P+A i — A+ P i ~ 
+ M 2 T R t^ 

+ 1+P ^ A / +lh EUx,Z : r>)u<,,X< 



(12) 



In the so-called chiral-odd sector in Eq. (fT2|) one may 
equally well work with T = <r' J+ = —e J T ia k+ ^. 

Notice that the GPDs for antiquaries are defined anal- 
ogously to the quark GPDs. One merely has to re- 
place in ([5]) the quark fields by the corresponding charge- 
conjugated fields. Unless stated otherwise all results in 
the following apply also to the case of distributions for 
antiquarks. 

As mentioned in the introduction we also want to con- 
sider possible relations for gluon distributions. The rel- 
evant correlation function for leading twist gluon GPDs 
reads 

F g M(x,A;\,\') 
1 f dz' 



xP^ 
x W ab 



2lT 



e* z <p';A'|F^(-H 
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z; ±z) F, 



+i{ 1 , 



|p;A) 



z+=0+ i 
zt — Ot 



(13) 



where the Wilson line is given in the adjoint representa- 
tion of the color SU(3), 



2 Z > 2 Z ) 



z + =0+ 
zt—Ot 



= [0+,-Iz-,0 T ;0+,iz-,0 T ] Qb 

= Vexp(-g dy-f abc At(Q + ,y-,OT))- (14) 

V J-(l/2)z- J 

The gluon field strength tensor in (fT3"|) has the standard 
form 

Fr(x)=d»A»(x)-d»AZ(x)+gf abc A b l (x)A»(x), (15) 

with f abc being the structure constants of the ££7(3). For 
the definition of the chiral-odd gluon GPDs we will need 
the symmetry operator S defined through 



SO ij = \{O lj + O jl - 5% o mm ) 



(16) 



for a generic tensor O lJ . One readily observes that the 
symmetrized tensor S O 13 has only two independent com- 
ponents. The twist-2 gluon GPDs are given through the 
correlator in Eq. (fT3"|) according to 

F 9 (x, A; A, A') = 6% F 9 ^ (x, A; A, A') 
1 



2P+ 



«(p',A') 7 + H 9 {x,t,t) 



l(J 



E g {x^,t) )u(p, A), 



2M 

F 9 (x,A;A,A') = it^ F 9 ^ (x, A; A, A') 
1 



(17) 



u{p',X')[j + l5 H g (x,^t) 



2P+ 

F*' ij (x,A;A,\') = -SF g M(x,A;\,\') 



(18) 



2P+ 



p+Ai _ A + pi 

T T -u{p',\')[ia +3 H 9 T {x^,t) 



+ 



2MP^ 

3 

T 



7 +A J 



A+ 7 ^ 



2M 



P+A 3 - A+P 3 - 

\p T H 9 T {x^t) 



+ 



7 +P£ - P+- 3 



— E 9 ( 



E 3 T (x,Z,t))u{p,\). (19) 



Note that the definitions of the chiral-even quark and 
gluon GPDs directly correspond to each other [compare 
the right-hand side (RHS) of (JTOJ) and lfIT|). as well as 
the RHS of (JTTJ) and fig])]. On the other hand, in the 
chiral-odd sector the definitions of the quark and gluon 
GPDs are (symbolically) connected by 



n 



A+P 



2MP+ 



— e^F^ k 



(20) 



[compare the RHS of (fT2)l and (HU)]. We also mention 
that our definition of all gluon GPDs X 9 differs by a 
factor x from the one of Ref. , 



= -X 9 {x^,t) 

here X 



Ref. 



(21) 



The main advantage of this choice in the context of our 
work is that the structure of the relations between GPDs 
and TMDs for quarks and gluons will look alike. 

Altogether there exist eight leading twist quark GPDs 
and eight leading twist gluon GPDs. All GPDs are real- 
valued which follows from time-reversal. Moreover, using 
commutation relations for the parton fields, they obey 
symmetry relations of the type 



X^ 9 (x,U)=±X^ 9 (-x,U), 



(22) 



where the minus holds for all GPDs X except H and E. 
It is therefore sufficient to consider only the region x > 



4 



as we will do in the present work. Eventually, hermiticity 
implies 



(23) 



where the plus holds for all GPDs X except Et- The 
relation (|23[) is needed later on in order to write down 
the general structure of the GPD correlator for £ = 0. 



B. GPDs in impact parameter space 

In a next step we want to consider the GPDs in trans- 
verse position (impact parameter) space pill in a. 

Of particular interest is the case £ = 0, where a density 
interpretation of GPDs in impact parameter space may 
be obtained [4l|. Such an interpretation, in principle, 
allows one to study a three-dimensional picture of the 
nucleon. In the context of the present work the impact 
parameter picture is relevant for different reasons: first, 
the intuitive picture for various transverse SSAs in semi- 
inclusive processes given in (35l . |36| is based on the im- 
pact parameter representation of the GPD E q . Second, 
the quantitative relation between the Sivers function /^? 
and the GPD E q , obtained in the framework of a scalar 
diquark spectator model of the nucleon |37l] , also con- 
tains E q in impact parameter space. Third, the impact 
parameter representation was used to point out analogies 
between chiral-odd quark GPDs and TMDs [Hj]. Fourth, 
we use this representation as a guidance to obtain new 
possible relations between GPDs and TMDs in the gluon 
sector. 

For the following discussion it is convenient to intro- 
duce a state describing the incoming nucleon with both 
longitudinal and transverse polarization as a superposi- 
tion of states with definite light-cone helicity (i 



\p; S) = cos(itf) \p; +) + sin(±0) e iv \p; -) 



(24) 



If one transforms this state to the rest frame of the nu- 
cleon, it describes a particle whose three-dimensional spin 
vector is given by 

S = (S 1 , S 2 , S 3 ) 
= (smd cos ip, sin$ smip, cos$) , (25) 

see Ref. (38|. In the following we use the notation 
A = S 3 = cos By means of the definition in (|24[) a cor- 
responding state (p'',S\ for the outgoing nucleon (with 

the same spin vector S) can be specified. Replacing now 
the helicity states in the correlators for the quark and 
gluon GPDs according to 



(p';X'\^(p';S\, \p;X)^\ P ;S) 



(26) 



the matrix elements for all possible helicity combinations 
can be obtained by an appropriate choice of the spin vec- 



tor. This statement is obvious looking at the relation 



F(x,A;S) 
= i[F(x,A 



+)+F(x,A,-,-)] 
±\[F(x,A,+,+)-F(x,A,-,-)] 
±S T [F(x,A,-,+) + F(x,A,+,-)} 
IS T [F(x,A,-,+)-F(x,A,+,-)] 



(27) 



which can be readily verified. Before considering the 
transformation to the impact parameter space we also 
give the definition of the light-cone helicity spinors. 
The calculations are performed using the conventions of 
Ref. [H, 



/ V2p 



u(p, 



^2V 2 p+ 



Pt ' 
V2p^ 
Pt ~ 



u(p,~) = 



V^ 7 



2 n + 



V 

/ -Pt- 
s/2p+ 

Pt ~ 
\-V2p^ 



ip\ 



ip\ 
V ip T 



ip\ 



\ 



(28) 



(29) 



When defining the GPDs in impact parameter space we 
restrict ourselves to the case £ = for which a density 
interpretation exists [4lJ. For convenience we also use 



Pt = 0t- These two conditions imply p + = p 



1+ 



P+ 



and A + = A =0. The parton correlators in impact 
parameter space are now given by the Fourier transform 



F(x,A T ;S). (30) 



The impact parameter 6t and the transverse part of the 
momentum transfer At are conjugate variables. In the 
impact parameter representation one naturally obtains 
diagonal matrix elements which is the crucial prerequi- 
site for a density interpretation. This gain of the impact 
parameter picture becomes evident after introducing the 
states (M El H 



|P+,6 T ;5) =Af 
(P+,b T ;S\ = TV* 



d 2 pT 

d 2 p' T 

(2tt) 2 



~^- %T \p-S), (31) 



(p';S\ 



(32) 



which characterize a nucleon with momentum P + at a 
transverse position 6t and a polarization specified by S. 
The normalization factor Af in these formulas is given by 



d 2 p T 
(2tt)2 



(33) 



and therefore infinite. However, using wave packets in- 
stead of plane wave states this infinity can be avoided jU 
E3- With the states in Eqs. ([31]) and the correlators 
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defining the GPDs of quarks and gluons can be rewritten 
as 



(P+,0 T ;S|<A(zi)r 



F q[ri (x,b T ;S) 

2 J 2ir 

x W(zi;z 2 ) ij;(za) |P+,0 T ;5), 



(34) 



1 



with 



x Wab(zi;z 2 )F+ l (z 2 ) \P + ,0 T ;S), 



(0+,t^-,&t) 



(35) 



(36) 



Zl/2 - V" i-r 2 ' 

Obviously, the two correlation functions in (j3~4"|) and ([33]) 
are diagonal. 

In analogy with Eq. ([30]) we define the GPDs in impact 
parameter space according to 



x( Xl b 2 ) 



(2nf 



X(x,0,-A 2 T ). (37) 



Using this definition one finds after straightforward alge- 
bra that the correlators in Eqs. (fTt?]) ~ (fl"2"|) for the quarks 
and in Eqs. (| 1 T[) — (fl~9|) for the gluons, written in impact 
parameter space at the kinematical point £ = 0, take the 
form 



M 



= \H q/9 (x,b T ). 



(38) 
(39) 



M 



£ T (x,b%) +2U q T {x,b^ t 



6| 



+ S T \H q T (x, 6|) - A b H q T (x, 6|) 
Ibip bp • St — Sp brp 



M 2 



H q T {x,b T ) 



(40) 



&>' ij (x,b T ;S) 



S b rp €.rj-i S rp 



M 



' ' ' l n T {x,b^)-^A b H T {x,b 2 ) 



S b l rpeip {^Ibp by • St — Sip bp 



M3 ^ ' ' 

In these equations we use the notation 



n 9 T {x^)\ .(4i) 



d 
db 2 



X(x,b%)), 



(42) 




FIG. 2: Kinematics for TMDs. 



P 



and analogous for the higher derivatives of the GPDs X, 
as well as 



A b X(x,b T ) = ^- 



b T 



db 2 



b 2 T ) 



(43) 



While Eqs. l|3"5)) - (|4"0"|) were already given in the litera- 
ture [H, Ull > the result in (|4"T|) is new. Since the point 
£ = is chosen, the GPDs E and Et do not show up in 
(J39J)— ([41]): the GPD E is multiplied by the kinematical 
factor A + = in the correlator, and Et vanishes due to 
the constraint in Eq. (|2"3"1) . 

The expression in (|38p , for instance, can be interpreted 
as the density of unpolarized quarks/gluons with mo- 
mentum fraction x at the transverse position bx in a 
(transversely polarized) proton. This density has a spin- 
independent part given by TL, and a spin-dependent part 
proportional to the derivative of £. Some details on the 
physical interpretation of (|39|) and (|40f can be found in 
Refs. BUI- 

Because of the spin-dependent term the impact param- 
eter distribution in (|38p is not axially symmetric (unless 
£' = 0), i.e., it depends on the direction of bx- In other 
words, the spin part causes a distortion of the distribu- 
tion P5]l . Note that the RHS in (|4"0"]) contains two terms 
providing a distortion, one determined by the first deriva- 
tive of Et + 2H.T and one given by the second derivative 
of Ht- In (HH) none of the three terms on the RHS is 
axially symmetric. Later on, we will use the results f55j) 
(|41[) and compare them with the corresponding correla- 
tors for TMDs. This procedure will give us some guid- 
ance in order to obtain possible relations between GPDs 
and TMDs (see also Ref. [38j]). It turns out that the 
specific form of the relations depends on the number of 
derivatives of the GPDs in Eqs. ([3"g ]) -(|H ]) . 



C. Transverse momentum dependent 
parton distributions (TMDs) 

In this subsection we summarize our notation for the 
TMDs. For the quark sector we follow the conventions of 
Refs. [5J, ly, l28|, 14711 . while for the gluons our treatment is 
similar to Ref. [4g|, where gluon TMDs were systemati- 
cally classified for the first time. In Fig. [2]the kinematics 
for TMDs is indicated. The nucleon momentum is de- 
noted by P (with P 2 = M 2 ), the quark momentum by 
k. The TMDs depend on k + , given by the longitudinal 



6 



momentum fraction 



(44) 



and on the transverse momentum kx- Like in the case of 
GPDs a possible polarization of the nucleon can conve- 
niently be described by a covariant spin vector S, whose 
components read 



S + = 



XP+ 

IT 



XP_ 

M ' 



(45) 



The quark TMDs are defined through the correlation 
function 



1 f dz" d 2 z T 
~ 2 J ~2^~ (2?r) 2 

x W +co (-±z;±z)</>(±z) |P;S) 



{P\S\$(-\z)T 



z+=0+ 



(46) 



in which a summation over the color of the quark fields 
is implicit, and only connected diagrams have to be con- 
sidered in pcrturbative calculations. Like in the case 
of GPDs, the dependence of the correlator in (j46|) on 
the renormalization scale is suppressed. The Wilson line 
in (|46[) is more complicated than the one for the light- 
cone correlators defining the GPDs [IE HH, E(| EE ■ 
It can be decomposed according to (see also Fig. [3]), 



W+oc(- 



2 Z > 2 Z ) 



z+=0+ 

\zt\ + , +oo~, - 
-±z T ;0 + ,+ocT 



i , 



", hz T ;0 



+ i . 



(47) 



where the future-pointing Wilson lines in (|47[) (running 
along the z~ direction in Fig. [3]) are appropriate for defin- 
ing TMDs in SIDIS In the Drell-Yanprocess the 
Wilson lines are necessarily past-pointing [21|, whereas 
in hadron-hadron scattering with hadronic final states 
even more complicated paths for the Wilson lines can 
arise [H, [H, EH, E|[ . Without loss of generality, in this 
work just the SIDIS case is considered. If different Wil- 
son lines in (|46[) are used, the results of our model cal- 
culations for the TMDs only differ by calculable factors. 

In general the use of lightlike Wilson lines in Ec|. (14"6")) 
can lead to so-called light-cone divergences [4!j EE l58| . 
and a refined definition of TMDs is needed in order to 
avoid this problem [23|, [2^, [4^, [5?], Ei| . This issue may 
make it more difficult to establish model-independent re- 
lations between GPDs and TMDs. On the other hand, 
for the nontrivial relations we are going to study in the 
context of our low order model calculations no light-cone 
divergence appears, and the use of the Wilson line in (147)) 
is safe. 



(Q+,\z-,\z T ) 



2 ' 2 
• «*- 



(0+, +oo-, \z T ) 
AM — t 



AM—* 

oo , -\zt) 



(^-\z--\z T ) (0 
FIG. 3: Path of Wilson line for TMDs in SIDIS 



The leading twist TMDs are obtained from the corre- 
lator in (|4"6"| by using the same three matrices T as in 

Eqs. CE0HE1, 



^(x,k T ;S) = ^ qh+] {x,k T ;S) 

— ff{ x > ^r) 
$> q (x,k T ;S) = ^ + ^(x,k T ;S) 



rp fVrpQrp , ^ 



(48) 



(49) 



T T h ±q/ 



■ h x q (x, k T ) + ——• h 1 ^(x, k T ) 



+ Srp ( h^rp^X, krp^j + 



2k 3 T k T ■ S T - S T k 2 . 



2AP 



(50) 



We note that a number of other notations exist for some 
of the quark TMDs (see, e.g., Refs. [1 El EE El), and 
that the corresponding TMDs for antiquarks are again 
given by the same correlation functions with charge- 
conjugated fields. 

The correlation function for the leading twist gluon 
TMDs reads 



§ 9[i3] {x,k T ;S) 
1 



xP^ 



dz_ d 2 z T k . z , | F+ j/_i, 
2n (2tt) 2 6 \ F > b \*« y 2' 



x W +00 , ab {-\z;\z) F+\\z) \P-S) 



2+=0+ 



(51) 



where the path of the Wilson line, now in the adjoint rep- 
resentation, corresponds to Eq. (j47|) . It is worthwhile to 
mention that (|5ip is not the most general gauge invariant 
operator definition for gluon TMDs. In the most general 
situation two (different) Wilson lines in the fundamental 
representation can appear (see also, e.g., Ref. [56]), and 
only in the particular case that these two lines coincide 
one obtains the definition in Eq. ([51) . However, corre- 
sponding to the discussion about more general paths for 
the Wilson line in the quark correlator (|4"6")l , also for the 
gluon TMDs a more general Wilson line in (|5"Tj) would 
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merely change the overall factor of our model results 
and not affect the general conclusions. For simplicity 
we therefore restrict ourselves to the definition (|5T|) . 

The twist- 2 gluon TMDs are now given through the 
correlator (|5Tj) according to 



$ 9 (x, k T ; S) = 5% 0, k T ; S) 



- fB 



fl ( x , k T ) 



trpl%rpjrp J_g 



M 



flT ( X 7^t) 



k T ; S) = iz% k T ; S) 



\g B 1L (x,k%) + 



Si 



A I 



9it( x > ^t) ■ 



(52) 



(53) 



$ T ' lJ (x, k T - S) = -S (x, fc T ; 5*) 

t-J fl/rp f\ rj 



2A/ 2 



h 1L (x, k T ) 



S fc T e J T fc (2fc T k T -S T - 5 T ±fl 



4M 3 



h$(x,h£). (54) 



Analogous to the GPD case, the definitions of the chiral- 
even quark and gluon TMDs correspond to each other 
[compare the RHS of (gHJ) and (JS2J), as well as the RHS of 
d49)l and (|53|) ]. In the chiral-odd sector the definitions of 
the quark and gluon TMDs are (symbolically) connected 
by 



T 



2M 



(55) 



[compare the RHS of $5$ and f5I])]. The gluon TMDs 
defined in Eqs. d52j) — ([54j) are related to those of Ref. (48] 



through 








/? = 


+G, 


hi 3 = 




-f^S _ 
J IT — 


— Gt , 


h 9 - 
IT — 


-[aH t 


a 9 - 
i/lL — 

i/lT — 


-AG L , 
-AG t , 


h ±9 - 
h ±9 - 

"it — 


-AHi, 



(56) 



In the forward limit ff corresponds to the unpolarizcd 
gluon distribution, and g\ L coincides with the gluon he- 
licity distribution often denoted by AG. Our notation 
for the gluon TMDs is a natural extension of the nomen- 
clature for quark TMDs. In the context of our work we 
prefer a notation with small letters for the TMDs in or- 
der to avoid confusion with the nomenclature for GPDs. 
Eventually, notice that yet another notation for gluon 
TMDs was proposed in Ref. [63j . 

Altogether there exist eight leading twist quark TMDs 
and eight leading twist gluon TMDs, which arc all real- 
valued. Using commutation relations for the parton fields 
one can derive the symmetry relations 



X , krp ) 



±Y q / 9 (- 



X , krp ^ 



(57) 



for the TMDs Y, where the plus holds for 
f U/a a/a h U h U h a h 

JIT ' WIT ' "l ' "lL ' "lT' " 



J-9 
IT > 



and the minus for 

fllB <ils 
J I 1 BIT ' 



u^a u±a h a h U 

"l ' "lL ' ''IT' IT 



(58) 



(59) 



Like in the GPD case it is therefore sufficient to consider 
only the region x > 0. 

The TMDs can also be divided into T-even and T-odd 
distributions. Whereas on the quark sector two func- 
tions are T-odd (/{^?, h^ q ), there exist four T-odd gluon 

TMDs h^l, h{ T , h$). 

The transverse momentum dependent transversity dis- 
tribution of quarks and gluons is given by 



i-Lq/g/ l2\ 

2 ll-^rp yX ; fcrp J 



2M 



(60) 



In the forward limit, where the parton correlators are 
integrated upon the transverse parton momentum, the 
gluon transversity drops out. This result is obvious from 
its prefactor in Eq. (|5"4")) . and ultimately follows from con- 
servation of angular momentum. 

With the exception of fx all TMDs characterize the 
strength of different spin-spin and spin-orbit correlations. 
The precise form of these correlations is given by the 
prefactors of the TMDs in Eqs. (|3g ]) -([50 j) and §2 § -$H $ . 
To be more specific, the TMD c/it and the transversity 
distribution in Eq. (f6T)|) describe the strength of a cor- 
relation between a longitudinal/transverse target polar- 
ization and a longitudinal (circular)/transverse (linear) 
parton polarization. By definition the spin-orbit correla- 
tions invoke the transverse parton momentum and either 
a polarization of the target (/fr), or of the parton (h^), 
or of both (giT, h^ Ll h^j). Note that a corresponding 
discussion on spin-spin and spin-orbit correlations also 
applies to the GPDs in impact parameter representation 
(see, e.g, Ref. [H[), where the role of the transverse par- 
ton momentum in the case of TMDs is played by the 
impact parameter. In this context the only difference 
between TMDs and GPDs lies in the fact that the spin- 
orbit correlations accompanied by the functions gix and 
h\ L have no counterpart on the GPD side for £ = be- 
cause the corresponding GPDs E and Et do not show up 
as explained after Eq. (|4"3"|) (see also the related discussion 
in Ref. [H ). 



III. RELATIONS BETWEEN GPDs AND TMDs: 
MODEL-INDEPENDENT CONSIDERATIONS 

In this section the current knowledge on attempts to 
establish model-independent nontrivial relations between 
GPDs and TMDs is summarized. As already mentioned 
earlier, a special role is played by the impact parameter 
representation of the GPDs, which was used in Ref. [36| 
with the aim to find a relation between the GPD E and 
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the Sivers function f^p. In [38[ the impact parameter pic- 
ture was exploited, in particular, to write down (model- 
independent) analogies between chiral-odd quark GPDs 
and TMDs. Here the treatment of Ref. [38[ is extended 
to the gluon sector. It is also argued that it is possible to 
consider GPDs in momentum instead of impact parame- 
ter space if one wants to have some guidance for possible 
relations between the two types of parton distributions. 
In addition, we briefly comment on the representation of 
parton distributions through light-cone wave functions, 
which also can give some hint on nontrivial relations be- 
tween GPDs and TMDs. 



B. Average transverse momentum 

In Ref. [36| an attempt was made to obtain a nontriv- 
ial relation between the GPD E q (in impact parameter 
representation) and the Sivers function fy^. The object 
considered there is the average transverse momentum of a 
quark in a transversely polarized target. The main steps 
of the treatment of [36| are repeated here. While in [36[ 
the light-cone gauge A + = was used, we do not work 
in a specific gauge. 

The average transverse momentum of an unpolarized 
quark in a transversely polarized target is defined by 



A. Forward parton distributions 

For completeness we first recall the well-known triv- 
ial relations between GPDs and TMDs. These relations 
hold due to the connection between GPDs and TMDs 
on the one hand and ordinary forward parton distribu- 
tions (PDFs) on the other. For the specific kinematics 
£ = t = some GPDs are related to the three twist-2 
PDFs: fi{x) (unpolarized distribution), gi(x) (helicity 
distribution), and h\{x) (transversity distribution) [60|. 
The same applies to some TMDs if they are integrated 
upon the transverse parton momentum. To be specific, 
in our notation one has 



f q/9 (x) = I d 2 k T f! /9 {x,k^) 



g\ /9 {x) = J d 2 k T9 l{ a (x,k3) 
= H q / g (x,0,0) = { d 2 b T n q/9 (x,b}), 



(61) 



(62) 



h\{x) = / d 2 k T ( h\ T {x, k 2 ) + ^ h±£(x, k 2 ) 



H q {x,0,0) 



d 2 b T (n q T (x,b 2 )-j^A b n q T (x,b 2 )). (63) 



Concerning the relations between TMDs and PDFs a 
word of caution is in order. Already in Sec. |TT]C we men- 
tioned that in general a refined definition for TMDs, con- 
taining certain non- lightlike Wilson lines, is needed to 
avoid possible problems with light-cone divergences. For 
such a definition the relations between PDFs and TMDs 
in Eqs. (|6ip -116 3 p. in g eneral, turn out to be nontrivial 
(see, e.g., Refs. |23ll57l.l59j). However, provided that one 
avoids the endpoints x — and x = 1, for our low or- 
der model calculations the definition of TMDs as given 
in (|46p and (l4"T)) does not lead to any problem and all 
relations (|6"Tj) - ([6"3")) are satisfied. 



(k¥(x)) UT = j d 2 k T ki r ^(x,k T -S) 
= 1J d 2 k T k* T [&(x,kT;ST)-$ q (x,k T ;-S T )] ,(64) 



with the correlator $ 9 in (|48j) . The second step on the 
RHS of (|64|) is justified because obviously only that part 
of $ 9 which depends on the transverse spin gives rise to 
a nonvanishing transverse momentum. In the following 
(kj^ '{x)) UT will be expressed in terms of (twist-3) quark- 
gluon-quark correlation functions [H, HH, [64| . To do so, 
we first rewrite the second term on the RHS of |64|) by 
means of the parity and time-reversal transformation, 

3> q (x,k T ;-S T ) 



2 J 2tt (2tt) 2 
x W +QO {-\z;\z)ip{\z) \P;-§t) 



1 f dz d 2 z T 



Ak-z 



■*(p-S T \i>(-\z)^ 



2 J 2ir (2tt) 2 
x W_oo(-±2;±z) 1>{\z) \P;S T ) 

This leads to the intermediate result 



z +=0+ 



(65) 



(k q /(x)) 



UT 



d k T k l T 



dz d 2 ZT 
~2T (2tt) 2 



x W +ao (-\z;\z) ~W-oo(-\z]\z) 
xil>{\z)\P;§r) 



.-+=0+ 



H7 + 



X W+oo(- \Z\\Z) 

XiP(tz)\P;S T ) 



zt—0t 



(66) 



In order to perform the second step in (|66[) we have ex- 
pressed the factor k l T through the transverse derivative 
d % T acting on the exponential, and performed an integra- 
tion by parts. 
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The derivative in (|66[) can now act either on the quark 
fields or on the Wilson lines. In the first case though, 
one gets no contribution to the average transverse mo- 
mentum, since the involved combination of Wilson lines 
vanishes, 



W +00 {-\z]±z) - W-oo(-±z; \z) 



,+=,,+ = 0. (67) 

zt—0t 



The result (J67J) is obvious because both Wilson lines are 
just running along the light-cone. 

On the other hand, if the derivative acts on the Wilson 
lines, one finds 



i 8 



v 2 Z ' 2 Z ) 



W+oo( 

= gf dy-W(-\z;y) t a F+ l (y)W(y;\z) 
= 2W(-±z;±z)I«> i (±z) 



z+=0+ 
zt—0t 



j,+ =z+=0+ 

Vt—zt—Ot 



zt—Ot 



(68) 



where the paths of the remaining Wilson lines run along 
the light-cone and the function I q ' L is defined by 



dy-\V{\z ]y ) t a F+\y)W{y;\z) 



y ' =z 1 
Vt=zt 



(69) 



Plugging the results together one arrives at the following 
expression for the average transverse momentum, 

(fcf; (x)) UT 

1 f dz 



,+=o+-(70) 

zt—Ot 



i J i r ^(P;5 T |^(-iz)7 + 
xW{-\ Z] lz)l^{\z)^z)\P;S T ) 



Equation (|7D|) is a representation of the average trans- 
verse momentum in terms of a sp ecific quark-gluon-quark 
light-cone correlator [36j, |52j, |64| . Since the gluon field in 
the three-parton correlator in (|70[) has zero longitudinal 
momentum one often talks about a soft gluon matrix el- 
ement. The reader is referred to [H, [(Jg, [§3, [68| where 
such (or similar) matrix elements were first discussed in 
connection with transverse SSAs. 

To unravel a possible connection between the Sivers 
effect and the GPD E q , in Ref. [H the RHS of fH 
was transformed to the impact parameter space, where 
it takes the form 



{k q r\x)) 



UT 



= r / d 2 bn 



dz 



— e™ p ^(P+,0 T ;5|^(z 1 ) 7 H 



x W(z x ; z 2 ) I q ' l (z 2 ) i>(z 2 ) \P+, T ; S) , (71) 

with z 1 / 2 as given in Eq. (|36p . Comparing the expression 
in (I7ip with the correlator ([M]) for the quark GPDs in 



impact parameter space (for T = j + ) one realizes that 
the only difference is the additional factor I q ' 1 and an 
integration upon the impact parameter br [36|]. On the 
basis of this observation one may hope to find a relation 
of the type 



(k$ x (x)) UT = J d 2 k T ki r <S> q (x,kT]S) 

<Pl*T l > i (x,l>T)J*(x,l>T',S), (72) 



where, in rough terms, the function X q ' 1 incorporates the 
effect of the gluon field in the correlator on the RHS 
of ([70]) . We mention that in the second term on the RHS 
of (|72|) only the spin-dependent term of T q contributes. 
Expressed in terms of TMDs and GPDs Eq. (|T2"j) reads 



UT 



M 



J d% T T q -\x,b T ) — 



jk,j qk / \ ' 



M 



(73) 



Interestingly, the relation (|T3[) is indeed fulfilled in 
the context of perturbative low order model calcula- 
tions [23] (see also Sec. Hvj) . It also provides an intu- 
itive understanding of the origin of the Sivers transverse 
SSA [aUSi]. However, Eq. {73]) does not have the sta- 
tus of a general, model-independent result (see also, e.g., 
Ref. (69|). The crucial problem lies in the fact that, in 
general, the average transverse momentum (k^ l (x)) UT 
caused by the Sivers effect cannot be factorized into the 
function 2 q ' 1 (called lensing function in [36[) and the dis- 
tortion of the impact parameter distribution of quarks 
in a transversely polarized target which is determined by 

(s q y. 



C. Generalization of relations 

To get further insight into possible relations between 
GPDs and TMDs, which at least may hold in the context 
of model calculations, we now follow a procedure given 
in Ref. [38j . The equations defining the GPDs in impact 
parameter space [see Eqs. ([3"8" ]) -pi ]) ] on the one hand 
and the TMDs [see Eqs. Ip ]) -((50] l and ([52 ]) -([53] l ] on the 
other obviously have a corresponding structure if one in- 
terchanges the impact parameter bx and the transverse 
momentum fey. Comparing these equations one directly 
finds out which functions may be related. However, using 
this procedure one cannot extract the precise form of the 
relations. Note also that the two TMDs g\r and h^ L have 
no counterpart on the GPD side, as already pointed out 
in SecHIlC. In the following we, respectively, talk about 
relations of first, second, third, and fourth type, depend- 
ing on the number of derivatives of the involved GPDs 
in impact parameter space. In the case of quark distri- 
butions the results given in this subsection were already 
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presented in Ref. [38(. At this point one has to keep 
in mind that, apart from the trivial model-independent 
relations (relations of first type), all relations presented 
in this and the following subsection so far have only the 
status of analogies between functions which follow from 
obvious analogies in the structures of the GPD and TMD 
correlators. Quantitative relations will be discussed in 
Sec. HVl in connection with model calculations. 

First of all, one finds the following connections by 
means of the mentioned comparison, 



fl' 9 <-> H q/9 , 



n q/9 , 



t y ' 



(74) 



which simply correspond to the trivial relations discussed 
in Sec. |TTT]A. 

Relations of second type contain GPDs with one 
derivative, 



£ q/g 



hp 



(s q + 2H q T ) 



h 9 1T 



2M 5 



h ±9 
"lT 



) » -( 



2(H 9 T 



IP 



(75) 



where the first relation in (|75|) involving and the 
derivative of £ q corresponds to Eq. (|73j) . At this point it 
is also worthwhile to notice that the computation of the 
average transverse momentum of a transversely polar- 
ized quark in an unpolarized target, using the correlator 
in Eq. ([50)1 . can be carried out completely analogous to 
Sec. IIIII B above where the transverse momentum caused 
by the Sivers effect is considered. Doing so, one even- 
tually obtains an equation corresponding to (|73|). with 
the quark Boer-Mulders function h 1 q showing up on the 
TMD side, and the first derivative of the linear com- 



bination £j, 



2T-0p on the GPD side. On the basis of 



these considerations one, in particular, also expects the 
same lensing function l q ' 1 to appear in the analogue of 
Eq. ([73]) . This feature indeed emerges in the context of 
the model calculations presented in Sec.|lV] We note that 
a corresponding discussion also holds for both relations 
between gluon GPDs and TMDs in Eq. (|75j). 

Finally, one obtains two relations of third type, con- 
taining GPDs with a second derivative, 



IT 



i(n q T )". 



hf 9 



2(£ 9 T + 2H 9 }j 



(76) 



and on the gluon sector even one relation where a GPD 
enters with its third derivative, 



i ±9 
'it 



(77) 



We emphasize that the number of derivatives of the 
GPDs in the various relations in a first place merely dis- 
tinguishes between the different type of relations, even 
though in specific model calculations one may find re- 
lations with exactly the number of derivatives showing 



up in (|74p - (|77|) . For the relations of first, second, and 
third type in Eqs. |(74]) - ([7d| this works, e.g., in the con- 
text of simple spectator models as outlined in Sec. IIVI 
below. We cannot check this feature for the relation of 
fourth type (|77|) by our model calculations because the 
respective parton distributions vanish. 



D. Relations in momentum space 

The relations presented in the previous subsection were 
obtained by comparing the equations defining the GPDs 
in the impact parameter representation (for £ = 0) with 
those defining the TMDs. Here we argue that there is ac- 
tually no need to make the Fourier transform to the im- 
pact parameter space. Relations corresponding to (|74[) - 
(l77|) also emerge if one compares Eqs. (|48t - ([50)l and (|52|) - 
([54")) on the TMD side with the momentum space corre- 
lators for GPDs in (UHJ) and (JT71) - (HH]) evaluated at 
C = 0. 

At the particular kinematical point £ = the RHS of 
Eqs. (Uni)— (H3) and dTTJ) — (JTHI) can be simplified consider- 
ably. Using the spin vector S introduced in Sec. [TT]B one 
finds after straightforward algebra 

F q / 9 (x,A T ;S) 
= H q / 9 (x, 0,-A|) - E q / 9 (x, 0,-A|) 



2M 



F q ' 9 {x,A T] S) 

= XH q / 9 (x,0,-A 2 T ), 
F q >i(x : A T ;S) 

= - l ^ L (E q T (x, 0, -A^) + 2H q (x, 0, -A? T ; 



,(78) 
(79) 



S J T (H q (x,0,-A^) 



A 2 

^T trf/ ft A 2 



AM- 



H q (x,0,-A' T ] 



2A - K -- 3 -- S ^H q {xA -A^, 



AM 2 



(80) 



F 9 ^(x,A T ;S) 



<> A* A- 7 
a i\ T i\ T { 



[£ T (x,0,-A^) + 2i^(x,0,-A^) 



2M 



H 9 T {x,Q,-A 2 T ) 



J,# T M,-A|) 



IS A T e^ (2A T A T ■ S T - S£ A 



8Af 3 
x H|,(a;,0,-Ay). 



(81) 



We repeat that the GPDs E and Et do not show up in 
the expressions above because of the choice £ = 0. 

One readily observes that the structure of Eqs. (178]) 
(j8Tj) corresponds to the structure of dS])-® and 
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(|54| . To be more specific, if one replaces — jAt/2 in 
the prefactors of the former equations by the transverse 
momentum kx one recovers the various prefactors of the 
latter. This strong similarity of the correlators for GPDs 
(in momentum space) and TMDs also allows one, like in 
Sec. IIIII C. to write down relations between the two types 
of parton distributions. In the following we refer again 
to relations of first, second, third, and fourth type, where 
the relations of first type are given by 

(Kt + Sp h^) ~ (#f + S H$) ■ (82) 
The relations of second type read 

fir' 9 «-» - ( Eq/9 ) > h i q ~ - ( e t + 2 #t) . 

(h\ T + J* hiS) - "2(ff# + £f*H), (83) 

and one finds 

hiS <-> 2 (#t) > h i 9 <->2(e% + 2H%\ , (84) 

for the relations of third type, as well as 

hiS^-^(H 9 T ). (85) 

for the one relation of fourth type. One may wonder why 
we distinguish here between different types of relations. 
The reason is simply that, in the defining equations, the 
prefactors of the GPDs and TMDs carry different powers 
of At and br, respectively. If there appears no At or &t 
in the prefactor we talk about a relation of first type, if 
they appear linearly we talk about a relation of second 
type and so on. In connection with this discussion it may 
be worthwhile to mention that the number of derivatives 
of the GPDs in impact parameter representation matches 
with the power of At in the corresponding prefactors 
in Eqs. (|T5|) - (f8"Tj) . This correspondence just reflects a 
property of the Fourier transform from the momentum 
space to the impact parameter space. 

Observe now the close analogy between the relations 
in (fg2" ]) -([g5 1) and those in (f73| ) — (f77 j) . In the context of 
our model calculations in Sec. IIVI we will consider rela- 
tions between TMDs and GPDs in both impact parame- 
ter space as well as momentum space. 

Eventually, let us briefly comment on nontrivial rela- 
tions between GPDs and TMDs that can be obtained 
by looking at a light-cone wave function representation 
of these objects. This issue was, in particular, dis- 
cussed in connection with a possible relation between 
the quark Sivers function and the GPD E q (see, e.g., 
Refs. [H [H [TO, (nl)- One finds that for both parton 
distributions the same light-cone wave functions appear, 
which hints at a close relation between / 13 ? and E q . How- 
ever, it turns out that in the case of the Sivers function 
one has to augment the involved wave functions by a 



phase factor not present in the case of E q [70]. This 
phase spoils a model-independent one-to-one correspon- 
dence between fi T and E q . So far only in the context 
of simple models the Sivers function can be represented 
through real wave functions multiplied by some addi- 
tional factor [ID, |4(| , allowing one to establish a direct 
connection between the two distributions. 

Moreover, in general (in full QCD) the representation 
of parton distributions in terms of light-cone wave func- 
tions may be problematic. This feature was already ob- 
served in the case of the unpolarized parton distribution 
fl(x) in Ref. @ . 

IV. RELATIONS BETWEEN GPDs AND TMDs: 
RESULTS OF MODEL CALCULATIONS 

While the previous Sec. IIIII is devoted to model- 
independent considerations on possible relations between 
GPDs and TMDs, here we discuss the relations in the 
context of specific model results. To do so, two specta- 
tor models are studied: first, the scalar diquark spectator 
model of the nucleon (see, e.g., Ref. [22j]); second, a quark 
target model treated in perturbative QCD. Some details 
concerning these models as well as the full list of GPDs 
and TMDs, computed to lowest nontrivial order in per- 
turbation theory, can be found in the appendices. 



A. Relations of first type 

We start this discussion with the relations of first type 
in Eqs. ([5T jl -(|o3 ]) , As an example consider the TMD /f 
in (|A7| and the GPD H q in (|A18|) in the scalar diquark 
model. Using these results one readily verifies that 

ft( X ) = J d 2 k T m x ,kS) 

_ g 2 (l -x) r 2 £ fcj + (m q + xMf 
2(2^)3 J T [k*+M*{x)Y 

= H q (x, 0,0) = J d 2 b T H q (x,b}), (86) 

i.e., that the relation (fuTj) holds. The unpolarized distri- 
butions of scalar diquarks and of quarks and gluons in 
the quark target model satisfy Eq. (f6Tj) as well. More- 
over, the results in both models obey the relations |62|) 
and (|63|) . Of course the model results have to satisfy 
these model-independent relations. Therefore, the com- 
parison discussed here merely serves as a consistency 
check of the calculation. 



B. Relations of second type 

In order to study the relations of second type in 
Eq. (|75|) in the context of our model calculations we first 
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consider the average transverse momentum of a parton 
caused by the Sivers effect (see Sec. IIIII B V Details of the 
calculation are given for a quark in the scalar spectator 
model. It turns out that the transverse momentum can 
indeed be expressed according to Eq. (f73|) in terms of 
(£ q )' ■ The reader is referred to (3?| where the connec- 
tion in (|75|) in the context of the spectator model was 
presented for the first time. 

By definition the average transverse momentum of an 
unpolarized quark in a transversely polarized target is 
given by [see also ([73]) ] 



/ 

g 2 e q e s (l — x)(m q + xM) 



UT 

r JkiJ ok 

= d 2 k T k T ^^f^( X ,k 2 ) 



d 2 kn 



fc| [k 2 +M 2 (x)] V M 2 (x 
g 2 e q e s (l — x){m q + xM) f d 2 l T l T 



2(2tt) 3 



(2tt) 2 f| 



d 2 k T 



JkjJ qk 



[k 2 + M 2 (x)] [(kr + h) 2 + M 2 (x)] 



(87) 



where we have inserted the result for the Sivers function 
from Eq. (|A8|) and used in the last step that 



(2tt) 2 f2 [(k T + T T y + M 2 (x)] 



Air fc| 



In 



k$ + M 2 (x) 
M 2 {x) 



If one now replaces the integration variable kx according 
to kT —> — kr — It one gets 



(k q /(x)) 



UT 



- 2 e q e s (l - x)(m q + xM) f d 2 T T V 



4(2tt) 
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Jkjj qk 

Klrp LrpiJrp 



[k 2 + M 2 {x)} [(kr + h) 2 + M 2 (x)} 

d I'J 1 Irp £rp IrpSrp 



~ J (2tt) 2 f2 (1 - x)M 



E q 



(^0,- IT § F ),(89) 



where we took for E q the result in Eq. (|A19|1 . In the next 
step the transformation to the impact parameter space 



is performed which leads to 

\^T ( X ))[JT 

d 2 lT 1>t 
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(27T) 2 f| 



jkjj qk 

rj-i trji >J r~p 
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(2tt) S 
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= / d 2 b T Z 



T X SDM 



(x, b T ) 



M 



£ q (x,b T ) 



(90) 



In (|90p integration by parts is used in order to obtain a 
representation of the transverse quark momentum that 
contains the derivative of £ q . Equation (|9"0|) coincides 
with the relation in (|73p between the Sivers function and 
the GPD E q , where the so-called lensing function [36| 

^sdm is S iven b Y 



^SDmO^t) 
_ e q e s (1 - 



(,r, 



(2tt) 2 



. / T 7 - b-T 



4tt 



(91) 



The calculation in ([57| ~ (f9T)|) also goes through step by 
step if one considers the average transverse momentum of 
diquarks in the diquark spectator model, using the Sivers 
function f^ T s from (jAl~6|) and the GPD E s from (|A25|) . 
One finds that the lensing functions for quarks and di- 
quarks in that model are the same, 



-SDM 



(92) 



Moreover, the transverse momentum (k T (x)} UT can also 
be computed for quarks and gluons in the quark target 
model. On the basis of the respective results given in 
Appendix |B] one obtains in that case the lensing functions 



T QTM( X ' b T) 



_ t9. 1 
— QTM 



(x,bq 



3g 2 (1-x) b T 



8tt 



b 2 



(93) 



Therefore, all results in the two spectator models satisfy 
the relation between the Sivers function and the GPD E 
in (|75j) . where the specific form of the relation is given 
in Eq. ([731 . In the context of the model calculations the 
lensing function does not depend on the parton type but 
merely on the model, which means in other words that it 
just depends on the target. Note that in both models the 
lensing function has the same overall sign and is negative. 

At this point we would like to add some discussion on 
the intuitive picture of the physical origin of the Sivers 
effect presented in Refs. [35|, [36|. The starting point is 
the observation that the impact parameter distribution 
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of unpolarized quarks in a transversely polarized target is 
distorted [44j (see also the discussion in Sec.HIlB), where 
the strength of the distortion is determined by the deriva- 
tive of £ q . A suitable measure for the distortion effect is 
the flavor dipole moment defined by [33 ] 



d q ' 1 = / dx 



d 2 b T b\ 



F q (x,b T ;S) 



2M 



ij qj 



2M 



(94) 



with the correlator T q taken from Eq. ([38[) . The dipole 
moment in (|94p is determined by the contribution n q of 
the respective quark flavor to the anomalous magnetic 
moment of the target. Up to this point the considerations 
are model-independent. 

Using SU (2)-flavor symmetry and neglecting in a first 
approximation the contribution from all other partons to 
the anomalous magnetic moment of the proton and the 
neutron one obtains 



K u/p _ K d/r 



1.7, 



K d/p _ K u/n 



-2.0. 



(95) 



On the basis of these numbers one finds dipole moments 
of the order 0.2 fm for the light quark flavors in the nu- 
cleon. This value is quite significant in comparison to the 
size of the nucleon [351 ] . 

It is now natural to speculate that this large distortion 
should have an observable effect. In fact, in Refs. [351.1361] 
it was conjectured that the Sivers effect observed in semi- 
inclusive reactions is intimately related to the distortion 
of the quark distributions in a transversely polarized tar- 
get. Though quite plausible on a first look, such a con- 
nection in general is actually not obvious, because im- 
pact parameter dependent GPDs a priori do not appear 
in the QCD-factorization formulas of semi-inclusive pro- 
cesses. At the moment the only thing known for sure is 
the following: if one describes the Sivers effect (in, e.g., 
SIDIS) in the framework of a spectator model (to low- 
est nontrivial order), the GPD £ q enters because of its 
relation to the Sivers function discussed above. In other 
words, so far the relation between the distortion of the 
correlator T q in impact parameter space and the Sivers 
effect is only established in simple spectator models (see 
also Ref. |69j). 

In the spectator models studied in Ref. [37J and in the 
present work, the Sivers effect, quantified by the aver- 
age transverse momentum, factorizes into the distortion 
effect times the lensing function according to Eq. ([9"0")) . 
In SIDIS for instance the lensing function describes the 
influence of the final state interaction of the struck par- 
ton. It is worthwhile to notice that the lensing func- 
tions in Eqs. (f9"Tj) and ([9"3"| only depend on the variable 
&r = &t/(1 — x) representing the transverse distance be- 
tween the active parton and the spectator parton. In 
Ref. (3?| it was shown that the lensing function in (|9Tj) 
is exactly the net transverse momentum which the active 
quark acquires due to its Coulomb interaction with the 
spectator diquark, if it moves from the position (c , c 2 , 0) 



to the position (c 1 , c 2 , oo). Therefore, (k T ,r {x)) UT in (|90|) 
is nothing but this net transverse quark momentum con- 
voluted with the transverse position distribution ([3"5|) of 
unpolarized quarks inside a transversely polarized target. 

What the picture for the Sivers effect given in Refs. 35. 
US H3 predicts is twofold: first, the Sivers effect for 
up and down quarks in the nucleon should have oppo- 
site sign. To arrive at this conclusion one has to make 
use of the (model-dependent) factorized form of the av- 
erage transverse momentum in (|9"0"|) . and of the fact that 
the distortion of the distribution T q in Eq. (|38|) (in a 
model-independent way) is given by the anomalous mag- 
netic moment n q . The opposite sign then follows from 
the phenomenological numbers in (f9"5"|) . At this point it 
should be noticed that the different signs are also the 
outcome of an exact, model-independent analysis of the 
Sivers function in the limit of a large number of col- 
ors N c [73j . The second prediction concerns the abso- 
lute sign of the Sivers function. The negative sign of 
the lensing function (f9"Tj) . reflecting an attractive final 
state interaction of the struck quark, implies for instance 
fiT^ p < J36|. It is interesting that the existing analy- 
ses [29l l30l [3lK l32l 33. 1341 of the data on the Sivers effect 
in SIDIS [ll], [jjl, Us [HI are m accordance with both 
predictions. 

Before proceeding to the Boer-Mulders function for 
quarks we briefly address the gluon Sivers effect in the 
nucleon. One can show in a model-independent way that 
the gluon Sivers function is suppressed compared to the 
Sivers effect of quarks [29(. This result follows from the 
Burkardt sum rule for the Sivers function [74], [7j| , stat- 
ing that the average transverse momentum of unpolar- 
ized partons in a transversely polarized target vanishes 
when summing over all partons, and from the large-iV c 
result according to which the Sivers function of the light 
quarks is equal in magnitude but opposite in sign [73j. 
The suppression of the gluon Sivers effect has also been 
confirmed by recent phenomenological analyses [70l . |76| . 

Now we would like to discuss the relation of second 
type in Eq. ([75]) for the Boer-Mulders function h\ q . A 
measure of the Boer-Mulders effect is the average trans- 
verse momentum {k T ^(x)Y TU of a transversely polarized 
quark (with polarization along j-direction) in an unpo- 
larized nucleon, which is given by the correlator ^ji 3 in 
Eq. ([50")) . A calculation completely analogous to (|57|) - 
([90)) in the scalar diquark model yields 



TU 



X - J d 2 k T k T [$$*(x,br;S) + $$ s (x ) b r ;-S)] 



d 2 k T k T ^^hi q (x,k*) 



x ( £ T {x,b T ) + 2H q T (x,b T ) 



(96) 
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In order to perform the last step in ([96)1 we used the result 
for h ^ q fro m Eg. Kg ) and the result for E q T + 2ll T from 
Eqs. (|A22[) and (| A23|) . The lensing function showing up 
in coincides with the one in Eq. (|9"Tj) for the quark 
Sivers effect, which means that it is not only independent 
of the parton type but also of its polarization. This fea- 
ture is actually not surprising if one keeps in mind that 
the lensing function arises from the high-energy (eikonal- 
ized) final state interaction of the quark. Such an eikon- 
alized interaction is not sensitive to the polarization of 
the quark. The relation in Eq. (|96p also holds for the 
results of the respective parton distributions obtained in 
the quark target model, provided that one uses the lens- 
ing function in (f9"5|) . 



According to Eq. (|96[) the (model-dependent) physical 
picture of the Boer-Mulders effect is completely analo- 
gous to that of the Sivers effect discussed above [39] . The 
magnitude of the Boer- Mulders function is proportional 
to the distortion of the impact parameter distribution 
of transversely polarized quarks in an unpolarized tar- 
get in (|40p . This distortion is given by the derivative of 
£ T +2H T , where in the spirit of Eq. ([94]) the magnitude of 
the distortion can be quantified in a model-independent 
way through the object [39| 



,9 _ 



dz(£ T (z,0,0) + 2ff T (x,0,0)V (97) 



quark target model one finds 
(fcr ( x )Ylt 

= \ J d 2 k T k T Jfc (x, k T ; S T ) - ^ jk (x, k T ; -S T )} 



a jJ Jel ql 

d2kTk T 2M 



x ^hl T (x, k^) + ^it( x ' ^t) 

n 9 T (x,b T ) ~ ^ A b n T (x,b 2 ) 



(98) 



where the two TMDs are taken from Eqs. (|B17|) 
and ([B18]> . and the GPDs from Eqs. ([B28[) and ([B30| . 
We note that the last line in ([54")) (term proportional to 
h^) does not contribute to the average transverse mo- 
mentum in (|98p . which is a model- independent result. 
In Eq. ([9"B"]) the same lensing function as in (j9"3"|) appears 
again. Unlike the average transverse momentum caused 
by the Sivers or Boer-Mulders effect, at present nothing 
is known about the magnitude of the gluon transverse 
momentum in Eq. (|98|) . 

To summarize, all our spectator model results obey 
relations as indicated in (f75|) . The specific form of the 
relations is given in Eqs. (|9"0"1) . ([TO]) , and (f9"5| and is of 
exactly the same structure for all three relations between 
GPDs and TMDs in ([75]) . 



While the anomalous magnetic moment n q is known from 
experiment [in combination with reasonable model as- 
sumptions, see (195[) ]. so far no experimental information 
on k t is available. However, results on k t were obtained 
in the framework of a constituent quark model ffl\ and 
in lattice QCD [78j . Both studies indicate that in nature 
k t of the nucleon may be as large or perhaps even larger 
than n q , and that it has the same (positive) sign for the 
two light quarks. Analogous to the discussion for the 
Sivers function, these results for n q T . together with the 
attractive final state interaction of the quark encoded in 
the lensing function, imply that both h\_ u ^ v and h\ d ^ v are 
sizeable and negative [39|, H§| . We note in passing that 
model-independent large-A^ considerations also predict 
the same sign for the Boer-Mulders function of the light 
quarks in the nucleon 

Eventually, we consider the relation of second type 
in ([75]) for the chiral-odd gluon GPDs and TMDs. Anal- 
ogous to the above discussion for the Sivers and Boer- 
Mulders effect we now study the average transverse mo- 
mentum of a linearly polarized gluon in a transversely 
polarized target. By definition this object is given by 
the correlator <E> T ' 13 in Eq. . In the framework of the 



C. Relations in momentum space 

In the previous subsection an important role was 
played by the impact parameter representation of the 
GPDs. Now we want to present relations between TMDs 
and GPDs in momentum space, where we again focus on 
the relations of second type [see the relations indicated 
inEq. (83)]. 

For the following considerations it is convenient to in- 
troduce moments of the GPDs X, 

= / d 2 A T ( " 1 X(x, 0,-A) , (99) 

2M 2 J \2M 2 J v WW' v ' 

and moments of the TMDs Y, 

Y^(x) = J d 2 k T (jtyY{x,U). (100) 

The relations discussed below are relations between the 
moments of the GPDs and the moments Y^ of the 
TMDs. 
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On the basis of the diquark model results for /{j? 
in l|A8p and for E q in (|A19p one finds by explicit cal- 
culation 

f^ (n \x) 

g 2 e q e s {l-x) (m q + xM)M 2n - 2 {x) H_ n 



16(2tt) 2 2™M 2 ™- 1 sin(wr) 

e q e s H- n T(2 - 2n) 



2(2tt) 2 (1 - x) T 2 (l-n) 



E q{n \x). (101) 



The relation (jlOip holds for < n < 1, i.e., n is not 
necessarily an integer. In (jlOip the T-function enters as 
well as H n , the analytic continuation of the harmonic 
numbers for noninteger n given by 



U 1 Jo l-x 



(102) 



Observe the relative factor (1 — x) between the mo- 
ments of the Sivers function and those of E q in Eq. ()10ip . 
The same relative factor was also found for instance in 
Ref. (7lT | by a model-independent investigation on the 
basis of light-cone wave functions. 

Relation (|101|1 also holds for the respective functions 
in the quark target model, provided that one replaces the 
quark and diquark couplings according to 



(103) 



This difference between the two models corresponds ex- 
actly to the difference between the lensing functions in 
Eqs. (j9T|) and (pjl . 

Our results in both the scalar diquark model and the 
quark target model obey all relations given in (|55|) . The 
specific form of the relations between the different parton 
distributions coincides with Eq. (jlOip . 

In a next step (jlOip is evaluated for three particular 
values of n, where for the most interesting cases one ob- 
tains 
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J IT 
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(2tt) 3 (1 - x) 



-■q^s 



4(2tt) 2 (1 - x) 



E q(1 \x). 



(104) 
(105) 
(106) 



We would like to add some discussion on these results by 
starting with the relation in Eq. (|106|) . This representa- 
tion of the relation between the Sivers function and the 
GPD E q was in principle already given in Ref. (37[, and 
is equivalent to Eq. (|89|) of the present work. This can 
be readily seen if one keeps in mind that 



d% l l T l J T A(f|) 



2°T 



d z l T lT A{1't) , (107) 



lT2 



which holds for an arbitrary function A(lj,). The rela- 
tion (|106p looks simpler than the one in ([50)1 containing 



the GPD in impact parameter representation. On the 
other hand, Eq. (|90p provides an intuitive physical in- 
terpretation of the Sivers effect (see Ref. [37| and the 
discussion in Sec. ITVlB). 

The relation in Eq. (fP04"|) was already recently discov- 
ered in Ref. Starting from (jlOip one has to make 
use of the result 



r*(i-nr Eqin){x) = ~i Eq{xA0) (108) 



lim 



H- n T(2 



in order to arrive at the relation (|104p . The identity 
in (|108p is valid for the result of E q in the simple specta- 
tor models considered in the present work but may not 
hold in general. 

Eventually, the case n = ^ in (|105[) is briefly addressed. 
It is this particular moment which appears in a natural 
way in cross sections because E q is accompanied by a 
factor At in the correlator (flT))) . while the Sivers function 
in the correlator (|48|) is multiplied by kx- 

It is worthwhile to mention that all relations presented 
in this subsection, exactly like the relation in (f90|) . also 
allow one to fix (in a model-dependent way) the absolute 
sign of the Sivers effect of light quarks in the nucleon. 
(Compare the corresponding discussion in Sec. IIVI B.) 



D. Relations of third type 

The discussion in Sec. IIVI B and IIVI C was exclusively 
devoted to the relations of second type. Here we want 
to consider the relations of third type indicated in ([76]) 
and (|M|) . As already pointed out in Sec. IIIII C one can 
expect the specific form of possible relations of third type 
to be different from the relations of second type. This 
expectation is also supported by the fact that all the 
TMDs in and are T-odd, while in (JTU) and $M$) 
they are T-even. To the best of our knowledge an explicit 
form of a relation of third type does not exist in the 
literature. 

Analogous to the result in poip a general relation be- 
tween the moments of TMDs and the moments of GPDs 
can be established. To be specific the diquark spectator 
model results for in (|A14|) and for i? T in (j A23|) obey 



h^ n \x) 



5 2 (1 



nM 2 



\x) 



4(2?r) 
1 



2 n M 2n ~ 2 sin(mr) 
nT(4- 2n) ~ q {n) 

tlrr 



(2tt)(1 



r 2 (2- 



(109) 



which again holds for < n < 1. The corresponding 
quark distributions in the quark target model fulfill ex- 
actly the same relation. Notice that here no replacement 
of the type (| 103[) is needed when making the transition 
from the scalar diquark model to the quark target model 
because the T-even TMDs and the GPDs in Eq. (fTTOl) 
receive the first nonzero contribution at the same order 
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FIG. 4: (a) Lowest order diagram for T-odd TMDs in spectator model calculations containing the interaction of the active 
quark with the target remnant. The eikonal propagator arising from the Wilson line in the operator definition of TMDs is 
indicated by a double line. Note that only the imaginary part of the box diagram on the left-hand side (LHS) of the cut 
is relevant for the calculation of T-odd functions. The Hermitian conjugate diagram (h.c.) is not shown, (b) Lowest order 
diagram for GPDs in spectator model calculations. The topology of diagram (a) matches with the one of diagram (b) if the 
quark-spectator interaction, described by the lensing function I q ' 1 , is factored out. 



in perturbation theory. In addition, we find that those 
gluon distributions in the quark target model, which en- 
ter the relation of third type as indicated in (|84p . satisfy 
a relation with exactly the structure of (|109[) . 

As expected, the general structure of the relation 
in ([TUT]) and in (fTT)!)|) is different. Note that due to the 
Wilson line contribution to the T-odd TMDs, the prefac- 
tor on the RHS in (| 1 1 [) contains couplings which do not 
appear in (|109p . Moreover, the relative power of (1 — x) 
between the moments of the TMDs and of the GPDs 
differs for both types of relations. 

Evaluating (|109[) for three specific values of n one finds 



,-1-9 (o) 
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(i- x y 
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i-Lg(i) 

U 1T 
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(2tt) 2 (1-x) 2 



(.*) = 
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(27T) (1 - X) 



- 2 H^\x). 



(110) 

(111) 

(112) 



Equations (|110p - (|112p are the counterparts of the rela- 
tions of second type in (|104p - (|106p . 

Keeping in mind the discussion in Sec. lIIII C [see in par- 
ticular (|76p ] one may wonder if the relation of third type 
in (|109p can be rewritten such that the second deriva- 
tive of the impact parameter distribution Hj, shows up. 
This is indeed possible for arbitrary values of n. Instead 
of providing a general formula we limit this discussion 
to the particular case n = 1 in which the most compact 
and appealing result follows. To this end we exploit the 
model-independent identity 



6| 







(27T) (1 - X) 



,H^\x). 



(113) 



In (I113p integration by parts is used in order to perform 
the first step. Combining now Eqs. (|112p and (|113p one 
immediately obtains 



h^\x)^ I Skr^h^ixM) 



(114) 



Note that this relation has a strong similarity to the re- 
lations of first type in Eqs. ([GT]) (16"3")) . Exactly the same 
result holds for the relation of third type containing the 
gluon distributions [see ifTB"])]. i.e., 



K 1 >(x) = / d'k 



2M 2 



d ^ T 2M* 2 [ £ t^t) +m g T {x,bl)) . (115) 



E. Relation of fourth type 

Eventually, the relation of fourth type indicated in (1771) 
and (|85p is considered. In the framework of the quark 
target model calculation at lowest order such a relation 
is satisfied because both the TMD hyj, and the GPD 
iff, vanish [see Eqs. (|B18|) and (|B30p ], In order to ob- 
tain nonzero results for those distributions higher order 
diagrams have to be studied. At present one can say nei- 
ther if higher order results obey a relation of fourth type 
nor how the specific form of such a relation could look 
like. One can only speculate that a possible relation of 
fourth type may be similar to the relation of second type 
because in both cases a T-odd TMD enters. 



F. Higher order diagrams 

As already pointed out above so far nontrivial relations 
between GPDs and TMDs are only established if the 
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(a) (b) 

FIG. 5: (a) Particular higher order diagram for T-odd TMDs in spectator model calculations containing the interaction of the 
quark with the target remnants, (b) Topology of diagram (a) if the quark-spectator interaction is factored out. Diagram (b) 
cannot represent a Feynman graph for a GPD because there is a mismatch between the number of particles on the LHS and 
the RHS of the cut. 



respective parton distributions are computed in simple 
spectator models to lowest nontrivial order in perturba- 
tion theory. It is therefore natural to ask what happens to 
the relations if higher order diagrams in spectator mod- 
els are taken into account. Here we would like to briefly 
address this point by focusing on the relations of second 
type. 

Consider for instance the relation between the Sivers 
function and the GPD E as given in Eq. (SQ). [The fol- 
lowing reasoning equally well applies to the relations of 
second type in |J96]) and in ([95]) .] We recall that l[90]) 
represents a factorization of the Sivers effect into the dis- 
tortion of the impact parameter dependent distribution 
of unpolarized quarks in a transversely polarized target, 
given by the derivative of the GPD 8 q , times the final 
state interaction of the active quark, described by the 
lensing function X q,t . 

Pictorially this factorization of the Sivers effect is indi- 
cated in Fig. 2J In order to compute the Sivers function 
to lowest nontrivial order in a spectator model one has to 
evaluate the cut-diagram in Fig. [4ja) . If in this diagram 
the quark-spectator interaction (lensing function) is fac- 
tored out, the topology of the remainder coincides with 
the diagram in Fig. > which represents a lowest order 
Feynman graph for a GPD. Hence, it is at least plausible 
that a factorization as given in Eq. (|5U)) can exist. 

Suppose now that higher order diagrams are taken into 
account for the calculation of the Sivers function, where 
a particular graph is depicted in Fig. OJa). Factoring out 
the quark-spectator interaction in this diagram, one ends 
up with the topology shown in Fig. [5](b) . The diagram in 
Fig. H^b), however, does not represent a Feynman graph 
for a GPD, because the number of particles on the LHS 
and the RHS of the cut does not match. Therefore, one 
can expect that even in spectator models the relations 
of second type no longer hold as soon as higher order 
diagrams are considered. Of course, strictly speaking our 
qualitative discussion here only suggests a breakdown of 
these relations and does not provide a rigorous proof. In 
any case, in order to arrive at a definite answer on this 
question a field-theoretic definition of the lensing function 
is needed. So far no such definition exists. 



V. SUMMARY 

Over the last few years several articles app eared (most 
notably Refs. [H, M, H3, H M, M, ml) suggesting 
nontrivial relations between generalized and transverse 
momentum dependent parton distributions. The present 
work is dealing with this interesting topic and has mainly 
a twofold purpose: first, a review of the current knowl- 
edge on nontrivial relations between GPDs and TMDs is 
given. Second, the existing results on such relations are 
considerably extended. 

In the following the new results of our work are listed. 

(i) The correlator for the chiral-odd leading twist 
gluon GPDs in impact parameter space is written 
down for the first time [see Eq. (|41[) ] . 

(ii) A definition of leading twist gluon TMDs is pre- 
sented which is very similar to the one proposed 
in and completely analogous to the definition 
of quark TMDs in 0,11, [H, |||. 

(iii) In the spirit of Ref. [38] nontrivial rela- 
tions/analogies between gluon GPDs and TMDs 
are obtained by comparing the correlators for 
GPDs in impact parameter representation with the 
corresponding TMD correlators. This procedure 
allows one to distinguish between different types 
of relations/ analogies, however does not provide 
an explicit form of a possible relation. The type 
of the relation is determined, e.g., by the num- 
ber of derivatives of the GPDs appearing in the 
correlator in impact parameter representation (see 
Sec. |HT]C). In our terminology relations of second, 
third, and fourth type represent nontrivial connec- 
tions between GPDs and TMDs. For instance the 
relation between the Sivers effect and the GPD E 
proposed in [35], [H, H3] is called a relation of second 
type. 

(iv) It is argued that the momentum space represen- 
tation of GPDs is also suitable in order to find 
possible relations between GPDs and TMDs. (see 
Sec. HUD). 
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(v) The first calculation of all leading twist GPDs and 
TMDs in both the scalar diquark model of the nu- 
cleon and the quark target model is performed (to 
lowest nontrivial order in perturbation theory). 

(vi) All our model results, in particular, also those for 
the gluon distributions, obey the relation of second 
type [see Eqs. ((H , ®, and (9SJ)], where the spe- 
cific form of this relation was first given in Ref . [37| . 

(vii) New relations (of second type) for moments 
of TMDs and GPDs in momentum space (see 
Sec. |IV]C) are found, which contain results pre- 
sented in Refs. [13, EO] as limiting cases. 

(viii) The first explicit form for the relations of third type 
is given. All our model results fulfill this relation 
[see Sec.lIVID. in particular, the compact results in 
Eqs. (fTT4|) and iflTgji ]. 

(ix) The results for the gluon distributions in the quark 
target model (trivially) satisfy the relation of fourth 
type in the sense that all involved distributions van- 
ish at the order in perturbation theory considered 
in our work. 

(x) It is pointed out that the relations of second type 
are likely to break down in spectator models if 
higher order perturbative corrections are taken into 
consideration. 



In very brief terms the general status of nontrivial rela- 
tions between GPDs and TMDs can be summarized in 
the following way: so far no model-independent nontriv- 
ial relations exist and it seems even unlikely that they can 
ever be established. On the other hand, many relations 
exist in the framework of simple spectator models. The 
phenomenology and the predictive power of the spectator 
model relation between the Sivers effect and the GPD E 
works quite well. This is the only relation which so far 
has been tested to some extent. Additional input from 
both the experimental and theoretical side is required in 
order to further study all other relations between GPDs 
and TMDs. Future work will certainly shed more light 
on this interesting topic. 
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APPENDIX A: SCALAR DIQUARK MODEL 

This appendix contains some elements of the scalar diquark model of the nucleon (see, e.g., Ref. [22]) and, in 
particular, the results of the various parton distributions in that approach. Since in the diquark model not only 
quarks but also diquarks are considered as elementary fields we also provide the definition of GPDs and TMDs for 
scalar partons in a spin-i hadron. 

The Lagrangian of the diquark model reads 

£sdm(i) = *(a;) [i-fd^ - M) *(z) + i>(x) {i-fD qill - m q ) ${x) + <p*(x) {Df D s ^ - m 2 s ) <p(x) 

-\F^(x)F^(x)+g[$(x)i!F(x)<p*(x) + *(x)il>{x)<p(x)] , (Al) 

where "J denotes the fermionic target field, ip the quark field, and ip the scalar diquark field. The essential ingredient 
of the model is a three-point vertex between the target, quarks, and diquarks, with the coupling constant g. This 
framework allows one to carry out perturbative calculations. All the results for parton distributions given below 
contain the coupling g to the second power, which is the lowest nontrivial order. 

In its simplest form the diquark model is of Abelian nature, where both quarks and scalar diquarks couple to a 
photon field via the covariant derivatives 

D q i ^(x) = [d'* + ie q A»(x)]i>(x), D»<p(x)=[d fl +ie s Af i (x)](p(x), (A2) 

with the charges e q and e s , respectively. In this model the target has no (electromagnetic) charge reflecting the fact 
that in QCD a hadron is color neutral. This condition directly implies e q = —e s . The field strength tensor of the 
photon is defined in the standard way by 

(x) = d»A v (x) - d v A tl (x) . (A3) 
Eventually, we mention that the condition M < m q + m s has to be fulfilled in order to have a stable target state. 
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FIG. 6: Cut diagrams contributing to quark TMDs and GPDs in the scalar diquark model: (a) diagram for T-even TMDs and 
GPDs; (b) diagram for T-odd TMDs. Note that in the case of GPDs the kinematics of Fig. [1] has to be used. 
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FIG. 7: Cut diagrams contributing to diquark TMDs and GPDs in the scalar diquark model: (a) diagram for T-even TMDs 
and GPDs; (b) diagram for T-odd TMDs. Note that in the case of GPDs the kinematics of Fig. [T] has to be used. 



For scalar fields only two GPDs (H s , E s ) and two TMDs (f( , f^) exist. In analogy with Eqs. 
defines the leading twist GPDs for scalars according to 



and (fTUl) one 
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The correlator in impact parameter representation (at £ = 0) is given by 



T s (x,b T ;S) = H s {x,b%) + 



€.rp I)rp S rp 
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(A4) 



(A5) 



i.e., it coincides in its form with Eq 
$ s (x,k T ;S) 
= xP 



The definition of TMDs for scalar fields reads 



dz d 2 ZT 



ik-z 



2tt (2tt) 2 



■*(p ] S\v*{-\z)W +00 {-\z;\z)i P (\z) \P;S) 



i+=0+ 



fi(x, k T ) 



2\ t T h, T D T fa f ^ 



M 



(A6) 



and is analogous to Eqs. (|1d| and We also note that in the diquark model an obvious change in the definition 

of the Wilson lines in the different parton correlation functions appears: instead of the strong coupling g the charges 
— e q and — e s have to be used in the correlators of quark and scalar diquark distributions, respectively. 

The perturbative calculation of the TMDs and GPDs in the scalar diquark model is basically straightforward. The 
parton distributions are given by the diagrams in Figs. [H] and [7] In the following we just quote the final results 
without providing any details of the calculation. We begin with the TMDs, where our results for T-even functions 
are of 0(g 2 ), while the first nontrivial results for T-odd functions exist at 0(g 2 e q e s ). This difference appears because 
necessarily a contribution from the gauge link is required in order to generate a nonzero T-odd function. Note that 
to this order all TMDs for antiquarks vanish identically. For the eight TMDs of quarks and the two TMDs of scalar 
diquarks one finds 



5 2 (1 



xMf 



2(2tt) 3 [k 2 + M 2 (x)Y 
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In the above formulas we used the abbreviation 

M 2 (x) = (1 - x)m 2 + xm 2 - x(l - x)M 2 (A17) 



for some specific combination of mass terms. All results in (|A7j) - (|A16|) were already given in the literature. To be 
specific, the full list of T-even quark TMDs was computed in Ref. [79|]. The quark Sivers function was considered 
in [2l|, (see also and the Boer- Mulders function in [8l|, [Hj]. The result for the Sivers function of the scalar 

diquark can be found in Ref. [83,]. 

Now we proceed to the model results for the GPDs, where in all cases nonzero results are obtained at 0(g 2 ). Again, 
to this order all GPDs for antiquarks vanish identically. We limit ourselves to the case £ = which is sufficient for 
the purpose of our work. The results read 
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FIG. 8: Cut diagrams contributing to quark TMDs and GPDs in the quark target model: (a) diagrams for T-even TMDs and 
GPDs; (b) diagram for T-odd TMDs. Note that in the case of GPDs the kinematics of Fig. [1] has to be used. 
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D q SDM (x, A T ; k T ) = [(k T - |(1 - x)A T f + M 2 (x)] [(k T + |(1 - x)A T ) 2 + M 2 (x)} , 
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We refrain from performing the integration upon the transverse quark momentum in the results given in (|A18[) - (|A25[) 
because this step is actually not needed for studying the relations between GPDs and TMDs. Note also that in the 
case of H and H the kx integration leads to a (logarithmic) ultraviolet divergence, which is not a peculiarity of the 
diquark model but rather a well-known feature of light-cone correlation functions. For the other GPDs this integral 
is finite, which for E is directly obvious and in the remaining cases can be shown by explicit calculation. 



APPENDIX B: QUARK TARGET MODEL 



The second model we are using is a quark target, treated in perturbative QCD. Mainly for two reasons this approach 
is of interest in the context of our investigation. First, the non-Abelian three-gluon vertex now enters the computation 
of various parton distributions. Therefore, one has an additional check of relations between GPDs and TMDs, beyond 
what can be done in the Abelian scalar diquark model. Second, the quark target model allows one to investigate 
relations between gluon distributions. 

The Lagrangian of the quark target model is given by 



C Q tm{x) = ${x) (*yD M - m) i>{x) - \F^{x)F^ a {x) , 



(Bl) 



i.e., it coincides the QCD Lagrangian for the specific case of a single quark flavor. The covariant derivative in (|B 1 1) 
reads 

D» 4(x) = [0" - ig t a A» (x)} t/,(x) , (B2) 

while the field strength tensor is defined in Eq. (fT5|) . 

We now list the results for TMDs and GPDs in the quark target model. Like in the case of the diquark model of 
the nucleon here we also avoid giving any details of the calculation. The relevant diagrams for the quark distributions 
are depicted in Fig. [51 and those for the gluon distributions in Fig. [5] The T-odd TMDs receive the first nonzero 
contribution at 0(g i ), while the results for all other distributions (T-even TMDs as well as GPDs) are of 0{g 2 ). 

For the eight TMDs of quarks and the eight TMDs of gluons one finds (in the kinematic region < x < 1), 
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FIG. 9: Cut diagrams contributing to gluon TMDs and GPDs in the quark target model: (a) diagram for T-even TMDs and 
GPDs; (b) diagram for T-odd TMDs. Note that in the case of GPDs the kinematics of Fig.Q]has to be used, and that diagram 
(b) contains a Wilson line for gluons. The open circle on the upper end of the gluon lines indicates a special Feynman rule for 
the field strength tensor in the definition of parton distributions for gluons (see Refs. [49, 83]). 
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At the kinematical point x — 1 there exist additional contributions for some of the parton distributions above, which 
arise if the on-shell intermediate state in the cut diagram is just the vacuum. Those contributions are neglected here 
for simplicity. Note that three TMDs in (|B3p ~ (|B18|) vanish. It is likely that a higher order calculation provides a 
nonzero result for these objects. Some of the results in (|B3j) - (|B18| were already given in the literature. Treatments 
of the quark TMDs ff, g\ T , hfl, h\ T can for instance be found in [13, H^,[84, 85, 86]. The Sivers function for quarks 
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and gluons was computed in Ref. [83| . 

Eventually, we consider the GPDs in the quark target model. The results (for £ = and < x < 1) read 
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Calculations of the chiral-even GPDs for both quarks and gluons can be found in Refs. [871. 18a. [89l| . To our knowledge 
the results for the chiral-odd GPDs in ([B19 ]) - ([B30ll are given here for the first time. 
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